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Abstract: In the context of global change, it is essential to promote the rational development and 
utilization of land resources, improve the quality of regional ecological environment, and promote the 
harmonious development of human and nature for the regional sustainability. We identified land use/land 
cover types in northern China from 2001 to 2018 with ENVI images and ArcGIS software. 
Meteorological data were selected from 292 stations in northern China, the potential evapotranspiration 
was calculated with the Penman—Monteith formula, and reanalysis humidity and observed humidity data 
were obtained. The reanalysis minus observation (RMO, i.e., the difference between reanalysis humidity 
and observed humidity) can effectively characterize the impact of different land use/land cover types 
(forestland, grassland, cultivated land, construction land, water body and unused land) on surface humidity 
in northern China in the early 21% century. The results showed that from 2001 to 2018, the area of 
forestland expanded (increasing by approximately 1.80x104 km’), while that of unused land reduced 
(decreasing by approximately 5.15x10* km?), and the regional ecological environment was improved. 
Consequently, land surface in most areas of northern China tended to be wetter. The contributions of 
land use/land cover types to surface humidity changes were related to the quality of the regional 
ecological environment. The contributions of the six land use/land cover types to surface humidity were 
the highest in northeastern region of northern China, with a better ecological environment, and the lowest 
in northwestern region, with a fragile ecological environment. Surface humidity was closely related to the 
variation in regional vegetation coverage; when the regional vegetation coverage with positive (negative) 
contributions expanded (reduced), the land surface became wetter. The positive contributions of forestland 
and water body to surface humidity were the greatest. Unused land and construction land were associated 
with the most serious negative contributions to surface humidity. Affected by the regional distribution 
pattern of vegetation, surface humidity in different seasons decreased from east to west in northern 
China. The seasonal variation in surface humidity was closely related to the growth of vegetation: surface 
humidity was the highest in summer, followed by autumn and spring, and the lowest in winter. According 
to the results, surface humidity is expected to increase in northeastern region of northern China, decrease 
in northern region, and likely increase in northwestern region. 
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1 Introduction 


Land use/land cover change (LUCC) is the most direct signal of the impact of human activities on 
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the natural ecosystems on the Earth's surface. As an important link between human socioeconomic 
activities and natural ecological processes, LUCC has become a hot topic in the study of global 
climate and environmental change (Sterling et al., 2012; Mooney et al., 2013; Liu et al., 2014; 
Tuffour-Mills et al., 2020; Zhang et al., 2020; Sun et al., 2021). Studies have shown that LUCC is 
closely related to the circulation of surface materials and life processes and that at least one third of 
the Earth's surface is directly affected by human activities in the processes of exploiting land 
resources and implementing land cover change (Shao and Zeng, 2012; Tian et al., 2012; IPCC, 
2019). In addition, as a key process in terrestrial ecosystems, LUCC plays important roles in the 
material cycle and energy flow. It mainly affects biogeophysical processes such as surface albedo, 
soil fluxes and vegetation evapotranspiration. Additionally, land coverage is essential not only for 
the transfer of energy between the land and atmosphere but also in the maintenance and 
optimization of ecosystem services, and it is also regarded as an "indicator" in the study of global 
climate change (Braswell et al., 1997; Sun et al., 2000; Wang et al., 2012; Buitenwerf et al., 2015; 
Verbeeck and Kearsley, 2015; Wang et al., 2016; Fu et al., 2017; Gao et al., 2019). 

Scholars worldwide mainly use two methods to study the climatic impact of LUCC. One 
method involves scenario simulations, which simulate the climatic factors assumed in specific 
scenarios by constructing a model and then predict the possibility of a certain impact of LUCC on 
a global scale (Lai et al., 2016; Lin et al., 2020). For example, Hurtt et al. (2011) reconstructed the 
interannual variation characteristics of global LUCC through a global land-use model and 
explored the response of LUCC to the global carbon cycle and climate change. Many scholars 
have also assessed the impact of LUCC on global climate through multiple climate models or 
interdisciplinary models (Oleson et al., 2008; Pitman et al., 2009; Li et al., 2017; Wang et al., 
2021). The other method is sensitivity analysis, which mainly explores the regional climatic 
impact of LUCC in a local area over a short time through remote sensing image data integrated 
with sediment, fossil or tree ring data (Dumont et al., 2015; Fezzi et al., 2015; Weinert et al., 
2016). Studies based on this method focused mainly on precipitation and temperature. For 
example, many scholars (Yang et al., 2009; Cao et al., 2010; Bian et al., 2019) have shown that 
LUCC has a far-reaching impact on regional climate, especially on precipitation and temperature, 
mainly by affecting the surface energy budget, the Bowen ratio and the distribution of 
precipitation among evapotranspiration, infiltration and runoff. Under the background of global 
climate change, regional surface humidity has a direct and significant impact on vegetation 
coverage. If only natural factors are considered, the regional vegetation coverage is high, resulting 
in high surface humidity, and vice versa; this relation has been previously validated. However, due 
to the close and complex relationship between the spheres of the Earth, different vegetation cover 
types may also influence the characteristics of regional climate (e.g., evaporation, transpiration and 
the water retention of roots, stems and leaves). Consequently, different impacts on regional surface 
humidity may also be observed. Therefore, in terms of the study of regional surface dry-wet 
conditions for long time series, the influence of LUCC on surface humidity cannot be ignored. 

Northern China spans a vast area (4.76x10° km), accounting for approximately half of the 
whole country. The climate here shows obvious regional heterogeneity and transitional 
characteristics. Most of northern China is in the dry-wet transition zone, so the ecological 
environment is highly fragile. Therefore, northern China is sensitive to global climate change. On 
the basis of a comprehensive consideration of various climatic factors, we objectively explored 
the impact of LUCC on surface humidity and its regional heterogeneity in northern China to 
provide a theoretical basis for human beings to rationally develop and utilize land resources, 
protect the regional ecological environment, and promote the harmonious development of the 
regional social economy and ecological environment. 


2 Data and methods 
2.1 Study area 


Located north of the Qinling-Huaihe Line, northern China covers a vast area with the longitude 
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and latitude ranges of 73°33'—135°05'E and 35°00’—-53°33'N, respectively. According to regional 
differences, we divided the 15 provinces, autonomous regions and municipalities in northern 
China into three regions: northeastern region (Heilongjiang Province, Jilin Province and Liaoning 
Province), northern region (Inner Mongolia Autonomous Region, Beijing City, Tianjin City, 
Shanxi Province, Hebei Province, Shandong Province and Henan Province) and northwestern 
region (Shaanxi Province, Gansu Province, Ningxia Hui Autonomous Region, Qinghai Province 
and Xinjiang Uygur Autonomous Region). In addition, there are significant differences between 
regional dry and humid zones in climate and vegetation types, which can span four dry-wet zones 
(humid zone, semi-humid zone, semi-arid zone and arid zone) from east to west, and the 
corresponding vegetation types change from forestland to grassland and even desert. It can be 
said that northern China is a sensitive area of global climate change. 


2.2 Surface humidity data and processing 


The climate data used in this study were derived from the daily surface climate datasets provided 
by the China Meteorological Data Service Centre (http://data.cma.cn). By excluding sites with 
missing data, 292 meteorological stations with complete observation data from 1991 to 2018 in 
northern China were selected for a comparative study with land use data (Fig. 1). The main 
climatic elements included in this study were the daily average temperature, daily maximum 
temperature, daily minimum temperature, sunshine duration, 2-min mean wind speed, daily 
precipitation, average water vapour pressure and atmospheric relative humidity. For each of these 
variables, we interpolated the missing or abnormal data from a few sites using the mean value of 
the data from adjacent sites in the same period. Furthermore, the seasonal differences in surface 
humidity and natural vegetation distribution in the study area were compared and analyzed by 
dividing the time series into four seasons: spring (March, April and May), summer (June, July and 
August), autumn (September, October and November) and winter (December, January and 
February). It is worth noting that based on the latitudinal position of the Qinling-Huaihe Line and 
the feasibility of map operations, we defined the area north of 35°N as the northern China, i.e., the 
study area (Fig. 1), to reduce the error caused in statistical analyses based on administrative 
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Fig. 1 Location of the study area (northern China) and spatial distribution of meteorological stations. SAR, 
special administrative region. 
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divisions. Therefore, the meteorological data and land use/land cover data from the provinces, 
autonomous regions and municipalities considered in this study represent the changes in north of 
35°N in China. 

The surface wetness index (SWI) is an important indicator used to measure the dry-wet 
condition of regional surfaces; it characterizes the income and expenditure of water at different 
spatial and temporal scales according to the relationship between precipitation and 
evapotranspiration (Hulme et al., 1992; Luo et al., 2016; Han et al., 2019). Based on the existing 
research and considering the influence of regional climate conditions on the accuracy of the 
reanalysis method, we selected the Penman—Monteith formula that has the highest accuracy and 
clearest physical meaning in both arid and humid areas, to calculate potential evapotranspiration 
(Allen et al., 1989, 1998; Du et al., 2001; Mao et al., 2011). The Penman—Monteith formula is also 
the only standard method recommended by the Food and Agriculture Organization of the United 
Nations (FAO) for calculating potential evapotranspiration. The calculation process is as follows: 


0.408A(R, —G)+y 0 u,(e, —e,) 
j A + y(1 + 0.34u,) : 
SWI=P/ET), (2) 


where ET» is the potential evapotranspiration (mm/d); Rn is the net surface radiation of crops 
(MJ/(m*-d)); G is the soil heat flux (MJ/(m?-d)), which is approximately zero when calculating the 
daily potential evapotranspiration; y is a wet-dry constant obtained from a table (kPa/°C); T is the 
average temperature at a height of 2.0 m (°C); u2 is the wind speed at a height of 2.0 m (m/s); es is 
the saturated vapour pressure (kPa); ea is the actual water vapour pressure (kPa); A is the slope of 
the saturated vapour pressure curve (kPa/°C); SWI is the surface wetness index; and P is the 
precipitation (mm/d). The larger the SWI value is, the more humid the surface, and vice versa. 

In research on the impact of LUCC on climate, some scholars have proposed the "observation 
minus reanalysis" method (Kalnay and Cai, 2003). According to this method, we defined the 
atmospheric relative humidity (ARH) data observed at 1.5 m above the ground recommended by 
the National Meteorological Research Center as the observed humidity values and the SWI as the 
reanalysis humidity value. In the process of calculating the SWI, both a series of influencing 
factors of the ARH and the underlying surface were involved. Therefore, it is necessary to use the 
principle of the reanalysis minus observation (RMO) method to subtract the observed humidity 
(i.e., ARH) from the reanalysis humidity (i.e., SWI) in the study area; in this way, the difference 
of surface humidity (i.e., RMO) can be obtained, and the spatial heterogeneity of surface humidity 
caused by differences in the physical and chemical properties of the underlying surface can be 
fully considered. The RMO can effectively characterize the impact of different land use/land 
cover types on surface humidity. Therefore, we selected the RMO to characterize the surface 
humidity at 1.5 m above the ground. The calculation process is as follows: 

RMO = SWI- ARH, (3) 


where ARH is the atmospheric relative humidity observed at 1.5 m above the ground 
recommended by the National Meteorological Research Center (%) and RMO is the difference 
between the SWI and ARH (%), which is mainly used to characterize the surface humidity. 

Since the slope of the linear trend can characterize the change trend of different indicators, we 
fitted the RMO values of 292 meteorological stations in the study area during the period of 
1991-2018 using univariate linear regression, and obtained the change trend of the RMO. The 
corresponding formula used to calculate the change trend is as follows: 


n n 


Slope = = lel (4) 
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where slope is the slope degree of the RMO linear fitting; n is the time series; i is the year; and x; 
is the RMO value in year i. The value of slope higher than zero means that the RMO is increasing, 
while the value of slope lower than zero means that the RMO is decreasing. 


2.3 Land use/land cover data and processing 


The MCD12Q1 remote sensing image data used in this study were provided by the Google Earth 
Engine (GEE) platform with a spatial resolution of 500 m. The image data were preprocessed and 
interpreted by ENVI 5.3 and ArcGIS 10.2, and the main processing steps included band 
combination, reclassification, image clipping and confusion matrix accuracy verification (Jing et 
al., 2017). On the basis of the International Geosphere-Biosphere Programme (IGBP) 
classification, we divided the land use/land cover types in the study area into six categories: 
forestland, grassland, cultivated land, construction land, water body and unused land (bare surface 
such as desert and Gobi) by using the method of combining supervised classification with visual 
interpretation classification and referring to the land use classification method of the Chinese 
Academy of Sciences (Liu et al., 2010, 2014). Finally, based on the GEE platform, we used the 
confusion matrix method to evaluate the accuracy of land use/land cover data in 2001, 2010 and 
2018 (Dai et al., 2021). The overall classification accuracies of the land use/land cover data in 
2001, 2010 and 2018 were 98.54%, 98.68% and 97.83%, respectively, and the corresponding 
Kappa coefficients were 0.969, 0.972 and 0.955, respectively, which satisfy the accuracy 
requirements of this study. 

The land use/land cover transfer matrix represents an application of the Markov model in 
LUCC, and this method can comprehensively and specifically describe the quantity and dynamics 
of mutual transfer of different land use/land cover types in a region. This matrix is derived from a 
quantitative description of the system state and state transfer in the system analysis. In this paper, 
we calculated the transfer matrix of various land use/land cover types in each period through the 
ENVI supervised classification method to obtain the change and transfer area of different land 
use/land cover types. Then, we introduced the area growth index to explore the changes in 
different land use/land cover types and their influence on surface humidity in the study area. 
Finally, combined with the expansion and reduction of land use/land cover types in different 
regions, we made a simple judgement on the future trend of regional surface humidity: 


1 
A _ A! 
= t 


N _x100%, (5) 


where N- is the area growth index (%); A; is the area of the i® land use/land cover type in period 
t+1 (km?); and 4; is the area of the i" land use/land cover type in period t (km?) (Jing et al., 


2017). 


3 Results 
3.1 Relationship between LUCC and surface humidity 


3.1.1 Time series characteristics of surface humidity 


The temporal variations of the ARH and SWI in the study area from 1991 to 2018 are shown in 
Figure 2a. This comparison shows that the ARH in northern China declined slightly but 
maintained stable during the study period, and the ARH generally exhibited a slight drying trend 
in the study area. The SWI, however, showed obvious fluctuation characteristics; specifically, 
before 2000, the fluctuation trend of the SWI was mainly downwards, with a small fluctuation 
range, while after 2000, the fluctuation range of the SWI was larger and the upwards trend was 
obvious. Therefore, it could be concluded that the variation trend of the SWI was not consistent 
with that of the ARH, and the difference between them could be attributed to the interference of 
the underlying surface influencing factors on the SWI being much stronger than that on the ARH. 
Combined with Figure 2b, the variation trend of the RMO in the study area was highly 
consistent with that of the SWI. The RMO also showed a trend of first decreasing and then 
increasing during the period of 1991-2018, with the year 2000 as the inflection point. The 
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fluctuation trends of both were almost synchronous. Accordingly, from 1991 to 2018, the surface 
of northern China showed a trend of drying first and then wetting, and the degree of wetness in 
the later period was much higher than that in the earlier period. In general, under the background 
of the downwards trend of the ARH, the surface of northern China became wetter during the study 
period (1991-2018). After the beginning of the 21* century, regional surface humidification was 
particularly obvious, which may be related to the impact of LUCC. 
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Fig. 2 Moving average time series of the SWI and ARH (a) and RMO (b) from 1991 to 2018 in northern China. 
SWI, surface wetness index; ARH, atmospheric relative humidity; RMO, reanalysis minus observation. 


3.1.2 Contributions of different land use/land cover types to surface humidity 

The land use/land cover types exhibited significant regional differences in terms of the RMO 
values due to their unique physical and chemical properties. In the areas with poor land cover, the 
land surface was mostly exposed and the moisture in the topsoil was easy to lose due to the direct 
influence of wind and solar radiation and other factors. Thus, surface humidity in these areas was 
low. Conversely, in well-covered areas, most vegetation had a protective effect on the land 
surface. At the same time, due to the water retention and water holding characteristics of 
vegetation, surface water was not easy to lose directly. In addition, plant evapotranspiration could 
increase the local ARH and promote the regional water cycle to a certain extent, ultimately 
increasing surface humidity. To highlight the sensitivity of the difference among the RMO 
values to land use/land cover types, we adopted the divisional statistical tool to determine the 
contributions of different land use/land cover types to the regional RMO in northern China in 
the 21* century (Fig. 3). 
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Fig.3 Average RMO values for different land use/land cover types in different regions of northern China during 
the period of 2001-2018 


From the perspective of the entire study area, the contributions of different land use/land cover 
types to the RMO had remarkable differences in the same region. Forestland and water body had 
the most positive contributions to the RMO. However, unused land and construction land had the 
most negative contributions to the RMO. In view of the different partitions, due to the different 
geographical environments, the regional differences for the same land use/land cover type were 
exceedingly obvious. Among them, forestland in northeastern region had the largest positive 
contribution to the RMO (66.54%), followed by that in northern region (21.93%) and 
northwestern region (15.20%). Similarly, water body, which made a large positive contribution to 
the RMO, also showed the characteristic of decreasing gradually from northeast to northwest. 
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Nevertheless, the negative contribution of unused land to the RMO was the lowest in northeastern 
region (19.63%), followed by northern region (—24.81%) and northwestern region (—29.80%). 
Because of the strong influence of human activities, the negative contribution of construction land 
to the RMO was also large. In northern region and northwestern region, the negative contributions 
of construction land were —0.72% and —14.82%, respectively. It is worth noting that the RMO 
contribution from different land use/land cover types was the highest in northeastern region and 
the lowest in northwestern region. Except for unused land, the other land use/land cover types in 
northeastern region all featured positive contributions to the RMO. In northwestern region, except 
for water body and forestland, the other land use/land cover types displayed different degrees of 
negative contributions to the RMO. 

In summary, the contribution of the same land use/land cover type to the RMO varied under 
different ecological environments. For example, in northeastern region, where the ecological 
environment was better, the contribution of forestland was obviously greater than that in 
northwestern region, where the ecological environment was fragile. Similarly, the obvious 
negative contribution of unused land to the RMO was far lower in northeastern region than in 
northwestern region, indicating that the impact of LUCC on surface humidity was related to the 
regional ecological environmental quality. 


3.2 Impact of LUCC on surface humidity 


3.2.1 Spatial variation characteristics of land use/land cover types 

The combination of the spatial distribution of LUCC in northern China in the 21* century (Fig. 4) 
and the statistical areas of LUCC in the study area in the same period (Table 1) reveals that the 
shrinkage of grassland in northeastern region was serious, decreasing by approximately 1.04x10+ 
km? from 2001 to 2018. The main degradation areas were concentrated in northern Heilongjiang 
and southeastern Liaoning. The remaining land use/land cover types showed various degrees of 
expansion. Specifically, forestland expanded by 0.67x10* km? and was mainly concentrated in the 
eastern and central forestland of northeastern region. In northern region, the areas of forestland, 
cultivated land and construction land also expanded to various degrees, with increases of 
1.10104, 0.95x10* and 0.41104 km?, respectively. The areas of the other three remaining land 
use/land cover types decreased, with the area of unused land decreasing the most (1.53x104km?), 
followed by grassland (0.90x10* km/?). In addition, in northwestern region, except for the 
considerable reduction of unused land (decreasing by approximately 3.70x10* km7), the 
remaining land use/land cover types displayed different degrees of expansion, with the area of 
cultivated land expanding by approximately 2.23 10+ km? and the area of grassland expanding by 
approximately 1.36104 km. In conclusion, the ecological environment in northern China was 
improved in the study area, with the coverage of forestland expanding (increasing by 
approximately 1.80x10* km?) and the coverage of unused land reducing (decreasing by 
approximately 5.15x10* km?). 

3.2.2 Sensitivity of spatial variations in the RMO to LUCC 

The changes in the areas of different land use/land cover types across northern China in the 21* 
century were obtained using the transfer matrix (Table 2). The spatial distribution of LUCC (Fig. 
4) and the changes in the RMO trend values (Fig. 5) demonstrated a good correlation between 
them. In northeastern region, the areas of forestland and water body expanded with the largest 
positive contributions to the RMO, and the RMO trend was relatively large, with the most rapid 
growth rate of 0.34/10a. In northern region, especially in western Inner Mongolia and western 
Shanxi, where the areas of forestland and grassland expanded and the area of unused land 
reduced, the RMO showed an increasing trend, while the RMO in most other areas showed an 
overall decreasing trend because of the obvious negative contributions of unused land and 
construction land. Similarly, in northwestern region that was dominated by the reduction of 
unused land area accompanied by different degrees of expansion of forestland, grassland and 
water body, the RMO exhibited an increasing trend, and the surface of the region tended to 
become wetter. In some areas in which construction land and cultivated land expanded, such as 
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Gansu, Ningxia and Shanxi, the RMO trend value was the smallest (the minimum value of 
0.89/10a), and the drying trend of the land surface was the strongest. In general, in the 21* 
century, surface humidity in northern China mainly exhibited weak humidification. The change in 
surface humidity was closely related to the variation in regional vegetation coverage; when the 
regional vegetation coverage with positive (negative) contributions expanded (reduced), the land 
surface became wetter. 
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Fig. 4 Spatial distribution of land use/land cover change (LUCC) in northern China during the period of 
2001-2018. (a), expansion of land use /land cover types; (b), reduction of land use /land cover types. Blank areas 
in northern China are the areas with no change in land use /land cover types. 


Table 1 Area changes for different land use/land cover types in northern China during the period of 2001-2018 
Change in area (x10 km’) 


ee Forestland Grassland Construction land Cultivated land Water body Unused land 
Northeastern region 0.67 —1.04 0.14 0.09 0.06 0.08 
Northern region 1.10 —0.90 0.95 0.41 —0.02 -1.53 
Northwestern region 0.03 1.36 2.23 0.01 0.08 -3.70 


Note: The positive value indicates the area expansion, while the negative value indicates the area reduction. 


Table 2 Area growth index values for different land use/land cover types in different regions of northern China 
during the period of 2001-2018 


Area growth index (%) 


Region - a B e Cultivated Construction Water Unused 
es aac Forestland Grassland 
land land body land 

Heilongjiang 5.58 -9.41 1.86 6.96 18.70 69.47 

Northeastern region Jilin 1.73 -16.19 6.20 11.00 34.40 47.30 
Liaoning 2.73 12.98 -8.58 13.89 0.98 25.16 

Inner Mongolia 15.90 0.41 16.44 5.64 -12.45 —6.42 

Beijing and Tianjin -2.33 9.72 —12.64 13.13 4.89 21.42 

Northern region Hebei 20.53 2.77 —5.49 17.23 0.78 63.22 
Shandong 71.49 7.70 -2.82 20.87 1.16 38.05 

Shanxi 61.62 -15.89 35.63 6.17 39.67 16.29 

Shaanxi 2.39 -1.67 8.63 6.63 -16.55 -50.38 

Gansu 19.99 2.32 17.12 1.43 18.59 —4.83 

Northwestern region Ningxia 129.02 -2.58 21.69 17.85 -11.11 —61.37 
Qinghai 38.49 4.94 -13.17 1.52 9.89 -3.46 


Xinjiang —17.40 3.10 55.69 0.58 7.34 -3.09 
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Fig. 5 Spatial pattern of the RMO trend changes in northern China during the period of 2001-2018 


3.3 Sensitivity comparison of the seasonal RMO changes to land use/land cover types 


Northern China is located north of the Qinling-Huaihe Line where the temperature changes 
clearly in the four seasons (the average temperature being lower than 0°C in January) and the 
surface coverage varies significantly. Based on these characteristics, it is necessary to explore the 
variation in the RMO in different seasons in northern China. 

The spatial correlation between the RMO values in the four seasons (Fig. 6) and the natural 
vegetation distribution (Fig. 7) was analyzed and compared. The seasonal variation in the RMO 
values in northern China was closely related to the distribution of natural vegetation. In general, 
the RMO values decreased gradually from east to west in different seasons in the study area, 
which was in good correspondence with the pattern of more natural vegetation distributed in the 
east and less natural vegetation distributed in the west. Among them, the RMO values of the 
whole study area was low in spring, and the range of high values was not large, mainly distributed 
in northeast region where forestland is concentrated. In summer, the RMO values increased 
significantly, and the range of high values expanded from forestland in northeastern region to 
grassland in northern region and the eastern part of northwestern region. In autumn, the RMO 
values in the study area showed a decreasing trend, and the range of high values began to reduce 
compared with those in summer. In winter, the range of the high RMO values concentrated in 
forestland of northeastern region. At the same time, land surfaces of northern region and 
northwestern region entered an abnormal drought period. It is worth mentioning that in northern 
China, the low value centers of the RMO in the four seasons were always concentrated in unused 
land. 


4 Discussion 


The transformation of land use/land cover types is a dynamic process, and LUCC is the product of 
this process. With the rapid development of social economy, human beings are continuously 
reshaping the land surface of the Earth, and the area is getting bigger and bigger, which has 
changed the land use/land cover types and vegetation coverage to some extent. It is generally 
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known that China is one of the regions with the highest intensities of human activities in the 
world due to its complex terrain and diverse underlying surface types, therefore, it is also one of 
the most representative research areas of LUCC (Hua et al., 2014). Similarly, northern China not 
only spans a large area but also features complex and diverse terrain, so the east-west climate 
difference is obvious. Thus, affected by the uneven distributions of regional water and heat, the 
LUCC in each region presents the characteristics of synchronous expansion and reduction, and 
the ranges of change are quite different. 

LUCC is one of the most important ways that human activities affect the Earth's climate. 
According to the report of the National Research Council (2005), the contribution of LUCC to 
climate change will even exceed the contribution of greenhouse gases. In the end of the 20" 
century, Chinese governments have implemented ecological protection projects to improve the 
ecological environment, control soil and water loss and prevent wind and sand hazards in northern 
China (Zhang et al., 2009; Mao et al., 2017). Affected by these efforts, the land cover situation in 
northern China has changed greatly. On the one hand, the climate effect of global warming with 
increasing emissions of greenhouse gases has been widely noted. On the other hand, compared 
with humidity data, temperature data have the advantages of more extensive, easier to collect, 
longer time series, good continuity, etc. Therefore, studies on the impact of LUCC on the regional 
climate have mainly concentrated on temperature, whereas research on the regional degrees of 
dryness and wetness is relatively limited. Kaufmann et al. (2007) studied the impact of 
urbanization on precipitation in the Pearl River Delta and found that rapid urbanization has led to 
a reduction in regional precipitation. Similarly, studies have shown that since the 1980s, rapid 
urbanization in Beijing has been closely related to the decrease in local summer precipitation 
(e.g., Zhang et al., 2009). In addition, other studies suggest that urbanization has led to a decline 
in relative humidity in urban areas and a slight increase in relative humidity in suburban areas (Li 
et al., 2008; Zheng and Liu, 2008; Zhao et al., 2011). There is no doubt that the pattern of LUCC 
affects the local climate. Therefore, understanding the change trends of land use/land cover types 
allows us to not only assess human activities on land use types but also judge the direction in 
which land use/land cover types are changing to predict the development trend of regional surface 
humidity. It could be inferred from the transfer matrix (Table 2) that surface humidity in 
northeastern region will increase in the future, but it is necessary to pay attention to the protection 
of regional ecology and prevent the continuous expansion of unused land. Similarly, the RMO in 
northern region will decline in the future, and land surface here will tend to become drier; 
therefore, it should pay attention to control the expansion of construction land strictly. At the 
same time, in urban planning, the construction of artificial wetlands (artificial lakes, wetland 
parks, etc.) and forest parks and other humidification projects could be appropriately increased to 
improve the surface humidity and promote the suitability of the living environment to a certain 
extent. Finally, the RMO values in northwestern region may increase in the future, and land 
surface will become wetter. This region should continue to strengthen the governance of the vast 
desert areas and improve the quality of the regional ecological environment. 

Through the above analysis, it was found that the seasonal variation in the RMO values in 
northern China had a strong response to the growth of regional natural vegetation. In winter and 
spring, due to the limitation of climatic conditions, plant species in grassland could not grow 
normally, then the grasses dried out and the leaves withered. The original green vegetation 
basically lost the ability to retain water. Therefore, the physical and chemical properties of the 
underlying surface of grassland were similar with those of unused land, and surface water was 
easy to be lost directly. Therefore, in winter and spring, the land surface of most areas in northern 
China was dry except for that in forestland of northeastern region. The degree and range of 
surface drying in winter were much larger than those in spring. In contrast, the climatic factors in 
summer and autumn were very beneficial to the growth of regional natural vegetation. Especially 
in summer, the vegetation coverage in most areas of northern China reached the maximum for the 
year. On the one hand, the expansion of vegetation-covered area was conducive to the maximum 
absorption and retention of surface water by the vegetation, so surface water was not easy to be 
lost directly. On the other hand, the expansion of vegetation-covered area promoted the regional 
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hydrological cycle to a certain extent, so the regional ecological environment presented a 
constructive cycle and further played a role in humidifying the regional land surface. 

In conclusion, the ARH in northern China declined slightly but maintained stable during the 
study period, and the ARH generally exhibited a slight drying trend in the study area. This 
conclusion agrees with previous research results (Wang and Zhang, 2008; Bian et al., 2011; Yao et 
al., 2017). Under this climate background, we chose northern China as an example to analyze the 
change trends of regional land use/land cover types and surface humidity, and it was found that 
the feedback effect of land use/land cover types on surface humidity exhibits an obvious spatial 
heterogeneity. Finally, we assessed the change trends of surface humidity in different regions of 
northern China, and the conclusions are basically consistent with those of previous studies (Bian 
et al., 2013; Yuan et al., 2017; Jia and Zhang, 2019). However, by referring to previous research 
theories (Kalnay and Cai, 2003), we used the RMO method to simply obtain the surface humidity 
from the SWI; nevertheless, further analysis of the internal mechanisms and factors that strongly 
influence the two types of humidity is lacking. In the future, we should strengthen the positioning 
and research of the relevant data, carry out more numerical experiments and analysis, and 
distinguish the impact of LUCC on climate more accurately. 


5 Conclusions 


In recent years, the land surface in northern China has become wetter in general. The RMO values 
for different land use/land cover types featured significant regional differences due to their unique 
physical and chemical properties. In general, the change in surface humidity was closely related 
to the variation in regional vegetation coverage; when the regional vegetation coverage with 
positive (negative) contributions expanded (reduced), the land surface became wetter. In addition, 
the seasonal variation of regional surface humidity was closely related to the vegetation growth in 
different seasons under the influence of special climatic conditions in the study area. It can be 
inferred that surface humidity in both northeastern region and northwestern region of northern 
China will increase in the future, but the surface will tend to be drier in northern region. Thus, 
special attention should be paid to the rational planning of land use/land cover types to improve 
the quality of the regional ecological environment. On this basis, we will work on future studies 
to explore the planning status and problems of land use/land cover types in different regions and 
propose practical suggestions for the sustainable development of the regional ecological 
environment from the perspective of ecological restoration. 
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